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ABSTRACT 

Aims. In this paper we study the effect of the quadrupole-term in the gravitational potential of a rotationally deformed central 
Be star on one armed density waves in the circumstellar disc. The aim is to explain the observed long-term violet over red 
(V/R) intensity variations of the double peaked Balmer emission-lines, not only in cool Be star systems, but also in the hot 
systems like 7 Cas. 

Methods. We have carried out semi-analytic and numerical studies of low-frequency one armed global oscillations in near 
Keplerian discs around Be stars. In these we have investigated surface density profiles for the circumstellar disc which have 
inner narrow low surface density or gap regions, just interior to global maxima close to the rapidly rotating star, as well as the 
mode inner boundary conditions. 

Results. Our results indicate that it is not necessary to invoke extra forces such as caused by line absorption from the stellar 
flux in order to explain the long-term V/R variations in the discs around massive Be stars. When there exists a narrow gap 
between the star and its circumstellar disc, with the result that the radial velocity perturbation is non-zero at the inner disc 
boundary, we find oscillation (and V/R) periods in the observed range for plausible magnitudes for the rotational quadrupolc 
term. 
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1. Introduction 

Many Be stars exhibit long-term periodicities in the rela- 
tive intensity of the double peaked Balmer emission lines, 
the so-called V/R variations, which occur on tim escales 
of typ i cally a f e w yea rs to a few decades , see e.g. iDachsl 



l|l987|) . lHubeTt] lll994h a nd lOkazakil l)l997l) and references 
therein. Okazaki <ll99ll) has suggested that a one armed 
(m=l) oscillation mode in the circumstellar disc about 
the Be star is responsible for the V/R variations. In a 
near Keplerian disc a one-armed mode shows up as a 
very slowly revolving perturbation pattern: the locally en- 
hanced temperature and density causes stronger emission 
in one part of the disc. This results in either the violet- 
or red component of the double Balmer emission peaks to 
be brighter than normal when the compressed enhanced 
emission region is located in the approach ing, respectively 
receding part of the disc. lOkazakil lll99lft applied a point 
mass potential for the Be star and obtained a slowly pro- 
cessing global m=l pattern by the invoked pressure per- 
turbations. However, in this way one generates retrograde 
modes which tend to propagate throughout the disc, while 
observations showed the V/R variations to occur only near 



the central star. lPapaloizou et al 1 lll992l) improved on this 
by considering a natural extension of this model by includ- 
ing the quadrupole-term in the gravitational potential of 
the rapidly rotating and thus flattened Be star. Deviation 
from a X/r point mass potential is known to give rise to 
processing elliptic particle orbits. In a gaseous disc the 
deviation from a point mass potential leads naturally to 
slowly revolving prograde m=l oscillation modes (if ex- 
cited) which are co nfined to a region within a few stellar 
radii from the star. iPapaloizou et al.l Ijl992|) showed that 
the propagation region extends in fact from the inner disc 
boundary up to about the point where the modal pattern 
speed equals the local particle precession speed. The radial 
size of this natural oscillation cavity is of order the stel- 
lar radius R s as the quadrupole component decays rapidly 
with increasing distance from the star. 

Meanwhile observations have shown that the pertur- 
bation pattern in several Be star-disc systems with V/R 
variatio ns revolves indeed i n a prograde sense, e.g. in 
/3 1 Mon llTelting et al.lll994tl iHummel & Vranckenl l|l997|) 
and iHummeffelJaimscnTkr (fl997f) calculated line profiles 
from discs with revolving m=l modal patterns and showed 
that these eccentric modes can explain the observed vari- 
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ation of the asymmetric line profiles quite well. The ex- 
istence of prograde m=l modes close to the central star 
has also been directly detected by spectrally resolve d in- 
terfero metry of the Be st ars £ Tau l|Vakili et alJll998h and 
7 Cas l|Berio et al.lll999|) . 

lOkazakil l|l997h argued that, although the observed 
density patterns appear indeed to re volve in the pro- 
grade sense, the proposed mechanism bv lPanaloizou et alJ 
( f 992) does not work for the hot early type Be star sys- 
tems and he introduced a radiative force with the aim of 
confining the m=l density waves to the inner part of these 
hot discs and induce prograde motions. Okazaki is indeed 
corr ect that in calculations with a rigid inner disc bound- 
ary ijSavoniie k. Heemskerklll993j) unrealistically low disc 
temperatures h ad to be adopted. I n a m ore recent study 
of m=l modes iFift fc Harmanecl l|200fil) note sensitivity 
of model predictions to disc model parameters, but com- 
ment that the predictive power of the model is better for 
the cooler Be stars where radiative forces are not added. 
In the latter work power law disc density profiles and rigid 
boundaries were also assumed. 

However, in this paper we show, using both semi- 
analytic and numerical approaches, that if there is a nar- 
row low density region or gap between the star and its 
circumstellar disc, so that the disc no longer has a rigid 
inner boundary but one that is essentially freely moveable, 
and there is also a significant single surface density maxi- 
mum not far away, prograde m=l modes with periods of 
a few years can be obtained even for hot Be star systems. 
Although the mechanism by which a Be star ejects matter 
into the circumstellar disc is poorly understood, we note 
that for the mechanisms suggested in recent studies, such 
as non-radial oscillations in the rap idly rotating Be star 
l|Rivinius et aLlEflOll ISavoniiell2005l) . or magnetic effects, 
one may indeed expect such a gap region to form when 
material with specific angular momentum larger than that 
corresponding to equatorial break up rotation is impul- 
sively ejected from the surface during an outburst. 

The plan of the paper is as follows: in Sect. 2 we give 
the basic equations for the oscillation modes and consider 
them in the low frequency limit. In that limit we use varia- 
tional principles to show that rigid wall conditions are the 
least favourable for prograde modes and derive conditions 
on the stellar apsidal motion constant for their existence. 
These are found explicitly for a class of surface density 
profiles, with inner narrow low surface density regions ly- 
ing interior to nearby global surface density maxima that 
are also considered numerically. The semi-analytic and nu- 
merical approaches were found to be consistent. 

In Sect. 3 we go on to apply the model to a hot Be 
star: 7 Cas (BO. 5 IVe) and a cool Be star: 88 Her (B7 
Ve). In both these cases prograde modes with reasonable 
periods are found for an appropriate degree of rotatonal 
distortion and apsidal motion constant, without the need 
to introduce an extra radiation force. Finally in Sect. 4, 
we summarize our conclusions. 



2. Basic equations for m=l density waves 

Let us introduce a non rotating cylindrical coordinate sys- 
tem (r, p, z) with its centre coinciding with that of the Be 
star. We assume the circumstellar discs about Be stars to 
be geometrically thin (z << r), near Keplerian and also 
assume hydrostatic equilibrium in the z-direction (v z = 0). 
For simplicity we work with a two dimensional approxi- 
mation by integrating over the z-direction. The z-averaged 
equations of motion and the equation of continuity follow 
as 

dv r ^ dv r ^ Vcp dv r v 2 1 OIL d^f 
dt dr r dp r £ dr dr 



dv v 
~dt 



dv v 
dr 



v v dv v v r v v 



r dip 



1 an 1 a* 

£ r dip r dip 



<9£ + I d(Zrv r ) + 1 <9(£tv 
dt r dr r dp 



= (3) 



where we have introduced a vertically integrated sur- 
face density and pressure 



£ = 



pdz and II = 



Pdz 



(4) 



We adopt a polytropic (isothermal) equation of state 
II = K £ 7 with 7 = 1 and define a (constant) sound 
speed c = ^/dIT/d£ = \fK in the disc. 

Since the circumstellar discs in Be star systems have 
ncgilible mass, the gravitational potential W is that of the 
central Be star. It is essential to take the rotational de- 
formation of the rapidly rotating Be star into account. 
Assuming uniform rotation and only weak distortions a 
multipolc expansion for \fr may be limited to th e monopole 
and quadrupole terms l|Papaloizou et alJll992^ : 



*(r) 



GM S k 2 n 2 s 
r 3 r 3 



(5) 



where M s , R s and £! s are the mass, radius and angular 
rotation speed of the Be star and k 2 its apsidal motion 
constant l)Schwarzschildlll958]) . 

The unperturbed disc is assumed axisymmetric with 
negligible velocity in the radial direction, i.e. with velocity 
v(r) = (0, rO(r)), whereby the angular velocity f2(r) in 
the disc is given by the condition for radial equilibrium 



GNL 



r 



1 dn 
s d7 



(6) 



where we have used expansion (JSJ for the gravitational 
field and / = f2 s /fi c with Vl 2 c = GM S /R? S as a measure 
for the stellar rotation speed. 

2.1. The linearised equations for the density waves 

After substituting the Eulerian perturbations with har- 
monic time dependence 
t'(r,(p,t) = £'(r) e i( - at - m f\ 

v' r (r,ip,t)=v' T (r)e i (° t - m rt ; v' (r, <p, t) - v' (r) ^ , -») | 
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which describe a density wave with m-fold azimuthal sym- 
metry, in equations and neglecting all orders higher 
than linear we find, after applying the operator relations 
- (i m) , the linearized oscillation equa- 



tions 



dW 
"dr~ 



, 2mS] „. 
/! + W 



(8) 



dv' fl 1 dS m dh\ . . / m 2 a . „. 
dr L = -U + E dr" + F^d7 U^-c^) W(9) 



where we have eliminated v' v and introduced 
iy(r) = c 2 E'/E and a = a - mQ,{r) 



(10) 



the oscillation frequency in a locally corotating frame. The 
square of the epicyclic frequency is k 2 — (2il/r)(dh/dr), 
with h = r 2 tt. The angular rate at which an orbit would 
precess in the absence of a perturbation to the pressure 
force is given by 



Q — K. 



(11) 



From now on we adopt m — 1 to obtain one-armed density 
waves. 

2.2. Low frequency limit 

Although equations (JS] - |5J can be readily solved numeri- 
cally as they stand, we remark that the m = 1 modes we 
consider are low frequency modes satisfying \a\ <C f2. This 
situation arises because the modes correspond to ellipti- 
cal distortions that precess on a long time scale compared 
to the orbital period. In the low frequency lim it the equa- 
tions m ay b e reduced to a simp lified form as in lPanaloizoul 
()2002f) and IPapaloizoul l)2005(l . This enables general de- 
ductions concerning the periods of oscillation and the way 
they are affected by boundary conditions to be made. To 
adapt equations JS] - EJ to this limit we assume that the 
system is of finite extent and such that it can be char- 
acterised by a single dynamical time scale fig that is 
everywhere significantly shorter than the precession time 
i n „A nnn ,ttvU^ n — 0(CIq). We set 

is a rel- 



scale uj p 1 and such that we can write f2 



uj p = eu;p r , where e is a small parameter and u pr 
ative precession frequency scaled to be O{0,q). We define 
a new eigenvalue A also of order Q by setting a = eX. In 
addition we set W = ew and under the assumption that 
the disc is thin c 2 = ec 2 ,. 

Equations (JSJ and 10 then give to lowest order in e 



dw 
d7 

dv(. 
~d7 



-2i (A - w pr ) Vj 



2w 



1 

27 



1 d£\ , wfi 
w r + i-5-. 



E dr 



(12) 



(13) 



Here we have used the fact that to within a correction of 
order e, we have h = r 2 Vt = \/GM s r. We can eliminate w 
from equations ifT^ and ljT3"|l to obtain a single governing 



second order differential equation for Q 
form 



Y.r 1 / 2 v' r in the 



d 
dr 



(y.3/2 ^,2 

[_ M3 

EO 



dQ 
dr 



2r 3 / 2 (A 



pr 



)Q 



or in terms of the original unsealed coordinates 



d 
dr 



„3/2 c 2 

EO 



dQ 
dr 



2r 3/2 



LO p )Q 



(14) 



(15) 



Equation Ijl5(l and its associated eigenvalue problem is of 
Sturm-Liouville form. This enables us to make use of the 
associated v ariational principle for th e eigenfunctions and 
eigenvalues l|Courant k, HilberH ll953l) . We comment that 
it is possible to formulate this principle using the prod- 
uct of Q and an arbitrary function of r as perturbation 
variable rather than Q. However, the results will be the 
same. 

2.3. Variational Principle 

In order to formulate a variational principle for the eigen- 
values, we must first consider the boundary conditions. We 
consider the disc to have an inner boundary at r = ro and 
an outer boundary at r = r\. We also suppose that the 
surface density is small but non zero at the boundaries as 
necessitated by our adoption of an isothermal equation of 
state. Then we may consider general boundary conditions 
for H15|l of the form 

dQ 
dr 

At r = ro we take A = Ao and at r = r\ we take A — A\. 
At a rigid boundary, Q — 0, which can be thought of as 
corresponding to A — > oo with dQ/dr remaining finite. 
The boundary condition that the Lagrangian variation of 
the pressure vanishes, as used in the numerical work be- 
low, corresponds to A — (l/E)(dE/<ir). 

Recalling that for Sturm-Liouville problems the eigen- 
values are real and the eigenfunctions may be assumed 
to be real without loss of generality (Courant & Hilbert 
1953), we derive an integral expression for the eigenvalue a 
from equation l|15l) by multiplying it by Q and integrating 
over (rg , r\ ) in the form 



AQ 



(16) 



2SO 



r 3 ' 2 c 2 A 2 
2SO ^ 



r n r 3 / 2 Q 2 



dr 



where 



r 3/2 



lo p Q 2 



2£ft V dr 



dQ' 



dr 



(17) 



(18) 



and the subscript 0/1 denotes evaluation at the in- 
ner/outer boundary respectively. Note that when the 
eigenfunction decays rapidly with distance the outer 
boundary terms can become negligible and make no ef- 
fective contribution. This is in fact the case for the eigen- 
functions we later consider. However, for completeness we 
shall retain these terms in the general analysis. 
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We also recall from equations © and l(TT|) that correct 
to first order in e, 



r 2 dS . 



and to adequate accuracy Q may be replaced by the 
Keplerian value SIk = (GM^/r 3 ) 1 / 2 in the above. 

The expression (|17f) forms the basis of the variational 
principle for eigenvalues a. In particular, for finite Aq and 
Ax, the largest eigenvalue, associated with an eigenmode 
with no zeros, is an absolute maximum with respect to 
variation with respect to arbitrary functions Q. When a 
rigid wall condition is applied at a boundary, the same 
variational formalism applies, but trial functions nust be 
restricted to those that vanish at that boundary. 

In this paper we are concerned with the question of 
the existence of prograde modes which have a > 0. This 
issue is seen to be equivalent to that of whether the max- 
imum of (|17|l is positive with respect to evaluation with 
all permissible trial functions. 

2.4. The existence of prograde modes and the 
dependence on boundary conditions 

We here point out that the existence of prograde modes is 
boundary condition sensitive by firstly pointing out that 
there exists sets of boundary conditions for which the exis- 
tence of prograde modes is assured. We go on to prove that 
a rigid wall boundary condition applied at cither bound- 
ary minimizes the maximum of (jl7(l and thus is the worst 
form of boundary condition for enabling prograde modes. 

Consider the equation derived from (|15J) when ki = 
and (7 = 0. This is, with the help of (|19|) . found to be 



are available when a rigid wall condition is not used. In 
the latter case Q is unrestricted at the boundary. Thus, 
other things being equal, changing from a non rigid wall 
(19) con dition to a rigid wall condition can only reduce the 
maximum eigenvalue o. 

2.5. A condition on &2 for the existence of a prograde 
mode 

From the variational principle it follows that it is sufficient 
for a prograde mode to exist that the numerator of l|17|) 
be positive for some Q which leads to 



d 
dr 



dQ 
dr 



Qrd_ 
lTdr~ 



r 2 dE 
E~ ~ik 



(20) 



In principle a solution of i|2U|). which we denote by Q = 
Qo, can be used to define values of Aq and Ai to use 
for the boundary conditions. In fact A\ can be specified 
and then Aq determined by inward integration from the 
outer boundary. We comment that a specific example of 
a solution for a power law E cx r~ 7 is Qo oc r~ s , where 
(5=1+7/2+ ^l+ 7 2 /4. 

The boundary conditions, or equivalently values of Aq 
and A\ so determined, are such that when the general 
eigenvalue problem is solved, there will be a prograde 
mode with a > 0. This follows immediately from the fact 
that if the trial function Q = Qo is adopted, only the 
term proportional to ki survives in (|17|) which must then 
be positive. Hence the maximum eigenvalue a must also 
be positive. 

The fact that a rigid wall boundary condition always 
reduces a maximum of (|17|) follows from the fact that 
adopting a rigid wall boundary condition through restrict- 
ing trial functions to be such that Q = at one of the 
boundaries is equivalent to restricting the search for a 
maximum to a subset of the possible trial functions that 



r 3/2 



ujpQ 2 



r 3 / 2 c 2 A 
2Eft 




(21) 



Using (|19fl this leads to 



n -3/2, 



Q' 2 



dr > 




Q 2 dr 



r 3/2 c 2 A 



2EO 



rV2 c 2 A 

2Eft ' 



(22) 



This in turn leads to a variational principle for the mini- 
mum value of ki for which a prograde mode exists. This 
is just the minimum that can be found from l|22|l for any 
permissible Q. Again because the possible trial functions 
for rigid wall conditions are a subset of those allowed for 
other boundary conditions, the minimum is larger in that 
case as is the corresponding value of ki. Accordingly free 
boundary conditions are always more favourable for pro- 
grade modes than are rigid wall conditions. 

From general scaling arguments, we would expect that 
for surface densities sufficiently peaked near r — R s , that 
this leads to a sufficient condition of the form 



k 2 f > 



Cc 2 R s 
GAL 



(23) 



where £ is a constant depending on the details of the sur- 
face density profile and boundary conditions but which 
is larger for rigid wall conditions relative to for example 
free boundary conditions. A condition of the same form as 
(|23|l was found in the case of a simple example presented 
below. 

The variational principle for the minimum k2 for which 
prograde modes are possible is equivalent to the eigenvalue 
problem obtained from equation (|15(l when a is set to zero 
and /c2 is taken to be the eigenvalue. The minimum eigen- 
value gives the minimum value of &2 for which prograde 
modes occur. 



J.C.B. Papaloizou and G.J. Savonije: One-armed oscillations in Be star discs 



2.6. Density profiles that are highly peaked near the 
inner edge 

We briefly consider the application of (1221) to systems with 
surface densities that are highly peaked near the inner 
edge. To proceed we note that the trial function Q = £ 
satisfies the free boundary condition (|16|) . Using this and 
adopting a surface density that is such that it and its radial 
derivative are very small at the boundaries so that the 
boundary terms become negligible, after some integrations 
by parts, (12211 gives that it is sufficient for prograde modes 
to exist that 

k2 jy r ^f^)\r>_£^r. (24, 
or after using Kepler's law 

£V^( fe /^-|0*>O. (25) 

From the above it is apparent that if £ is highly peaked, 
with a maximum at r — R m ax say, that prograde modes 
are guaranteed if k 2 f 2 > c 2 R* nax /{GMR 2 s ). Thus for 
Rmax = l-6i? s , we obtain prograde modes if k 2 f 2 > 
4.1(c 2 R s ) / (GM) . We now consider a specific example of a 
family of profiles and evaluate the condition exactly. This 
turns out to be of similar form. Thus, as has also been 
confirmed by normal mode calculations, the general form 
of our results does not depend on details of the surface 
density profile. 



5 

Table 1. The values of the quadrupole factor k 2 f 2 above 
which there is a prograde m = 1 mode in the disc deter- 
mined from the semi-analytic theory as a function of the 
disc model parameters a and (3 defined in (|26|) . and the 
type of inner boundary condition. 



Inner disc boundary 


Q 





k 2 fGM s /(c^R 3 ) 


Free 


2 


7.5 


4.7 


Free 


2 


6 


3.9 


Free 


2.5 


5.5 


3.1 


Free 


2.5 


7.5 


4.0 


Rigid 


2.5 





17.3 



where r = 3\/4 + a 2 , and q = (6+3a+r)/18. Specification 
of an inner boundary condition, here taken at r — R s or 
2=1, then determines IC 2 and hence k 2 f 2 as an eigen- 
value. The smallest eigenvalue for fixed a and (3 then leads 
to the smallest value of k% for which prograde modes can 
exist. For a rigid wall condition Q = and for vanishing 
Lagrangian change in pressure, 3zdQ/dz — Q(a — 3/3), at 
2 = 1. 

Once an inner boundary condition is adopted, it 
is a simple matter to obtain the smallest value of 
k2f 2 GM s /(c 2 R s ) for which prograde modes are possible. 
Results for various a and (3 and different inner boundary 
conditions are given in tabled These are fully consistent 
with the numerical work we have carried out (see below 
for further discussion). 



2. 7. A simple example 

We consider a surface density profile with a form very sim- 
ilar to that adopted below for our numerical work, being 
given by 



£(/•) = .4 ( ^ ) <-x P 



-(3Rl 



(26) 



where A, a, and (3 are constants. 

For this form of surface density, the equation for Q, 
when a = is 



d 2 Q dQ (Z + a 3(3 
dx 2 dx 



x x 1 \x" x° 

Here x = r/R s and JC 2 = k 2 f 2 GM s /{c 2 R s ) + 3/3. In terms 
of z — x~ 3 , we find 

9z 2 + ^ (3(1 - a)z + 9/3z 2 ) =Q(-a- 2K 2 z) (28) 

Solutions of gSJ can be found as power series expan- 
sions. We adopt the solution that vanishes most rapidly 
as z — > 0, or equivalently x — > 00 which can be 
written, to within an arbitrary normalization factor, in 
terms of a confluent hyperg eometric function in the form 



terms 01 a confluent nypcrgC' 
JWhittaker fc Watsonlll flfiaji 

i-Fi(2/C 2 /(9/?) + «; r/9 + 1; -0z), 



2.8. Strength of the quadrupole factor A^/ 2 

The structure and dynamical behaviour of the disc 
depends on the strength of the quadrupole term in 
equation (JSJ which is proportional to the factor k 2 f 2 . 
Unfortunately, the apsidal motion constant k 2 and the 
stellar rotation speed factor / = f2 s /fi c are poorly de- 
termined stellar parameters. For B-type main sequence 
stars stellar evolution calculations yiel d for fc? typical val- 
ues in the range 2.5 x 10~ 3 — 10 -2 l|Claret fc Gimene^ 
ll99ltlClaret)ll995ft . The rotation speeds are derived from 
line width measurements which includes the uncertain 
inclination factor sini and the effect of rotational grav- 
ity darkening which is thought to have lead to too 
low estimates of the ro tation speeds in earlier studies 
llPorter fc Riviniusl 12003^ due to underweighting of the 
cooler, most rapidly rotating equatorial region. Yet it 
generally remains unclear how close a particular Be star 
is to critical rotation, although in some c ases interfer- 
ometry may be helpful l(Tvcner et al J 12004) and suggest 
large rotational deformation, e.g . for the Be star a Eri 



l|Domiciano de Souza et all 120031). This has been ques- 
tioned bv IVinicius fc Levenhagenl ((20061) . who deduce a 



(29) 



lower ratio of equatorial to polar radius for this star and 
a r otation rate Sl/S l c ~ 0.8. 

ICranmeil ( 20051 ) made a new detailed statistical study 
of the threshold rotation rates for the formation of Be 
star discs. He comes to the conclusion that the early type 
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Be stars (07e-B2e) exhibit a spread of equatorial rotation 
rates between a lower limit around 0.4 — 0.5 Sl c and an 
upper limit near break up speed £l c . He finds the spread 
in rotation rates to narrow progressively when considering 
later spectral types, approaching critical rotation speeds 

for the AOe stars . 

iFremat et all l)2005|) studied the effects of fast rota- 
tion in early type stars and calculated synthetic spectra 
with classical plane-parallel non-LTE model atmospheres 
for rapidly rotating stars. They too find that gravity dark- 
ening can lead to serious underestimates of Tasini val- 
ues when considering the rotational broadening of spec- 
tral lines. By comparing with spectra of 130 Be stars they 
conclude that the average angular rotation speed of Be 
stars is most probably il s — 0.9 f2 c . From all this we infer 
that for the hot, early type Be stars the quadrupole fac- 
tor /C2/ 2 may lie somewhere between 4 x 10~ 4 and 10~ 2 , 
unfortunately quite a wide range. 



3. Numerical results 

3.1. Adopted disc structure 

For the (poorly known) density distribution in the equilib- 
rium disc we adopt the same form (|26[) as used in the semi- 
analytic analysis above, except for an additional constant 
B which is chosen so as to provide a floor to the surface 
density such that it does not decrease to below 5 x 10~ 4 of 
its maximum value (and inhibit a solution for isothermal 
discs): 



E(r) = — exp( — 



B 



(30) 



The exponential factor with the parameter is introduced 
to cut off the steeply increasing density profile near the 
star and to produce a gap (or low density zone) just out- 
side the stellar surface. The formation of a gap seems a 
rather natural assumption for circumstellar discs that are 
generated by impulsive mass ejection from the central star. 
It allows us to apply a free boundary condition at the in- 
ner disc edge. As indicated above, the eigenvalue solution 
is only sensitive to the integrated disc properties and inde- 
pendent of the precise shape of the density maximum, as 
long as it is close to the star. Therefore simulations with 
an exponential cut off for different values of (3 should be 
adequate to probe the dependence of the results to the 
disc's density profile near the star. 

From here on we work in units for which G = M s = 
R s = 1 and, as for the semi-analytic discussion, take r = 1 
for the inner disc boundary ro, neglecting the small width 
of the gap between star and disc (the stellar radii being 
not well known anyhow). Observationally derived (crude) 
power laws (r~ n ) for the density distributi on p(r) in the 
circu mstellar discs are typically n ~ 2 — 4 l|Waters et all 
Il987ft . Adopting n = 3 or 4 we find for a semi-Keplerian 
disc with disc semi-thickness (in z direction) H ~ c/VL 
that a ~ 3/2 or 5/2. 



The profile of the surface density S(r) for a = 2 and 
(3 = 7.5 is plotted in Fig.^ This disc extends quite far out 
from the star to r ~ 40. The deviations from Keplerian ro- 
tation are strongest near the inner disc boundary where in 
our units Q = 1.0103 (1.000) and the epicyclic frequency 
K = 0.9881. At r = 10 the deviation from Keplerian ro- 
tation is already small with f2 = 0.0308 (0.0316) and 
K = 0.0307. The numbers quoted in brackets refer to 
the angular speed in a purely Keplerian disc (for which 

n = k). 

We study the influence of the disc structure on the ob- 
tained oscillation periods by also considering more com- 
pact discs with slightly different values of a and j3. We 
assume the temperature in the isothermal discs is given 
by Td — 2/3 T e , where T e is the Be star's effective temper- 
ature, unless stated otherwise. 

3.2. Boundary conditions and numerical solution 

When doing numerical work we do not make a low fre- 
quency approximation as above but solve equations JSJ- 
|5J) directly. However, as above, and unlike for previ- 
ous calculations of m=l modes in circum stellar discs 
l)Savoniie fc Heemskerklll993T: IOkazakilll997l) we consider 
the existence of a small low surface density region, or 
gap, between the star and surrounding disc. Formation 
of a gap may be expected when the disc material was 
ejected from the star as a result of non-radial oscillations 
( Rivinius et alJEoOll) or magnetic levitation with specific 
angular momentum larger than Q s R 2 S . 

We apply free boundary conditions to both the inner 
and outer boundary of the disc and require the Lagrangian 
pressure perturbation <5n to vanish at both boundaries, or 



v r = -77 



W 

2 dlnE 



(31) 



unless specifically noted otherwise. These conditions cor- 
respond to those used in the semi-analytic discussion. 
We integrate equations I|8I9|1 using a fo urth order Runge- 
Kutta method with adaptive stepsize l)Press et al.lll98fih 
away from both disc boundaries, adopting an arbitrary 
boundary value W=l, to a conveniently chosen fitting 
point in the central region of the disc. We then search for 
that (real) oscillation frequency a for which the difference 
in the ratio of W/v' r determined from integrations starting 
on different sides of the fitting point becomes sufficiently 
small. The linear solution can then be scaled on each side 
so that W and v' T match separately. We comment that 
for the cases considered below, the location of the outer 
boundary was chosen such that further extension did not 
affect the results. 



3.3. A hot Be star: 7 Cas (B0.5 IVe) 

We now consider models for 7 Cas, the brightest Be star 
in the northern sky for which we adopt M s = 15M0, ra- 
dius R s — 8R Q and effective temperature T e = 25000 K. 
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iFift fc Harmanecl l|2006|) adopted the same radius, but a 
larger mass (18.3M ) and higher effective temperature 
(T e = 33000 K) for this star. We note that, according to 
the semi-analytic result in Sect 12. 51 these larger values re- 
quire the quadrupole term k 2 f 2 to be larger by only a fac- 
tor oc c 2 /M s , i.e. by a factor ~ 1.08 in order to obtain the 
same oscillation period. The observ ed V/R quasi-period s 
of 7 Cas fall in the range 4-7 years l|Doazan et alJll98"^) . 

We first consider the quite arbitrary semi-Keplerian 
circumstellar disc model described by l|3T)l) with a = 2 
and (3 — 7.5, for which the density profile in the inner 
disc is shown in Fig. ^ Note that if we adopt a rigid inner 
boundary to the disc (i.e. adopt (3 — 0, so that there is no 
gap) with v' r — there cannot exist prograde m=l modes 
in the hot disc model considered here as the pressure forces 
dominate over the gravitational quadrupole effect and we 
obtain retrograde solut ions that propagate all over the disc 
( Papaloizou et al.lll992t) . 

However, if we assume there is a (narrow) low-density 
region or gap between the star and the disc and consider 
a free inner disc boundary with a non- vanishing radial ve- 
locity of the gas, the situation changes drastically. Then 
the gas in the disc can oscillate in the radial direction 
right from the inner boundary where the gravitational 
quadrupole term is strongest, so that the latter term can 
dominate the dynamic behaviour. The numerical calcu- 
lations show that a free inner boundary leads to a very 
different outcome whereby we do find prograde m=l os- 
cillations in the observed period range, see Fig. consis- 
tent with our semi-analytic estimates. An observed V/R 
period of 7 years would for the disc model with a = 2 and 
(3 = 7.5 thus require a rotation parameter k 2 f 2 = 0.0065 
which is rather large, although in the allowed range, see 
Sect. 12.81 Fig. ^shows the shape of the eigenfunctions for 
a m=l mode with period P = 6.6 years. The shape of 
the eigenfunctions varies little with k 2 f 2 (or period). The 
shape of the density perturbation with a peak near r = 2.3 
is in fact rather similar to what iBerio et al.l l(l999h have 
inferred from their interferometric observations of 7 Cas: 
they deduce a prograde revolving m=l pattern centred 
near r ~ 2.5 from the stellar centre. 

For the extended disc shown in Fig. ^ the also ob- 
served V/R period of 4 years would require a value of the 
quadrupole factor k 2 f 2 ~ 0.008, closer to the upper limit 
of 0.01. However, the oscillation periods of the m=l modes 
studied here appear sensitive to the size and adopted den- 
sity profile of the disc. An alternative profile £(r) corre- 
sponding to a — 2.5 and (3 = 5.5 is plotted in Fig. |3J to- 
gether with the m=l eigenfunctions for k 2 f 2 — 0.005. The 
maximum for £ has shifted slightly more inward compared 
to the profile shown in Fig. ^ while the density vanishes 
already near r ~ 9. Note that in this case the disc extends 
to only r ~ 9, i.e. just beyond the propagation boundary 
of the m=l mode near r ~ 6.5. For this more compact 
(measured by the location of the surface density maxi- 
mum) disc the oscillation periods for a given quadrupole 
factor k 2 f 2 are significantly shorter and extend to smaller 
rotation values, as can be seen in Fig. For a compact 
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Fig. 1. A model of 7 Cas with M s = 15M©, R s = 8R© and 
T e = 25000K. The disc profile is defined by a = 1.5 and 
(3 = 7.5 for which the disc has a free inner boundary (with 
a narrow gap between disc and star) as well as a free outer 
boundary near r ~ 40. Plotted are the density profile E(r) 
of the unperturbed disc together with the eigenfunctions 
£' and v' r , in fact 3 (v' r ), for the m=l mode with P ~ 7 
yr. We only show the inner 10 stellar radii from the stellar 
centre. The vertical scale is arbitrary. 
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Fig. 2. The oscillation (or V/R) periods in years versus 
the quadrupole factor k 2 f 2 of (prograde) m=l modes for 
the extended disc in the 7 Cas model with a — 2 and 
(3 = 7.5shown in Fig. ^ together with the results for the 
compact disc with a = 2.5 and (3 = 5.5 shown in Fig. |5J 
The observed V/R quasi-periods for 7 Cas range from 4-7 
years. To see the effect of a steeper decline of S towards 
the inner disc boundary, we have also plotted the results 
for (3 = 6. 

(unevolved?) disc like this an m=l oscillation period of 4 
years requires a smaller quadrupole factor of k 2 f 2 ~ 0.005. 
The m=l oscillation periods for an intermediate disc with 
a = 2 and (3 = 6 are also plotted in Fig. 

One can speculate that just after mass ejection and 
re-building of the circumstellar disc it may be rather com- 
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pact and exhibits V/R quasi-periods in the lower part 
of the observed range while after some time, when an- 
gular momentum transport has taken place in the disc, 
it may get smeared out leading to longer periods of the 
m=l modes. This type of evolution would in Fig. ^corre- 
spond to a vertical shift from the compact disc with the 
plotted lower branch (/? = 5.5) to the extended disc with 
the higher branch ((3 = 7.5). That the two branches are 
not fully consistent in that they require slightly different 
quadrupole factors to explain the observed range in V/R 
periods should not surprise us in view of the arbitrariness 
of our disc models. It is interesting that the evolution from 
a 4 to 7 year V/R period occurred during the period 1970- 
1990 which coi ncided indeed with the build up of the disc 
around 7 Cas (|Doazan et al.lll987| . 

It is also of interest to compare these results with 
those obtained using the semi-analytic theory. Values of 
k*2f 2 GM s /(c 2 R s ) above which prograde modes exist for 
the values of a and f3 used in Fig. [5] are given in ta- 
ble ^ These can be converted to values of k 2 f 2 using 
GM s /(R s c 2 ) — 1.5 x 10 3 which is appropriate to this 
model. These values correspond to the asymptotic limits 
for P -> 00 in Fig.El They are found to be 0.0031, 0.0026 
and 0.0021 for the pairs (a = 2, @ = 7.5), (a = 2, (3 = 6) 
and (a = 2.5, (3 = 5.5), respectively, and are accordingly 
consistent with the plots in Fig. [3 
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Fig. 3. A model of 7 Cas with a compact disc defined by 
a = 2.5 and f3 = 5.5. The circumstellar disc extends only 
to r ~ 9, with a steeper decay of E towards the inner 
boundary in comparison with the disc shown in Fig. ^ 
Plotted are the density profile E(r) of the unperturbed 
disc with the eigenfunctions E' and v' r for the m=l mode 
with period P ~ 4 yr. The vertical scale is arbitrary. 



3.4. A cool Be star: 88 Her (B7 Ve) 

We now turn to models for 88 Her, a Be star with es- 
timated mass of M s — 5M Q with R s = 4R Q and 
T e = 12500 K. For this cooler Be star system we consider 



two different discs, one with a = 2.5 and (3 — 7.5 as before, 
i.e. with a gap between star and disc, whereby the disc ex- 
tends out to r ~ 19, and another model with /3 — 0, so 
that the inner disc touches the stellar surface. The bound- 
ary layer that exists in that case (assuming il s < Q c ) 
between the star and the inner disc is ignored here and 
replaced by the condition of a rigid wall: v' r = 0. The re- 
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Fig. 4. A model of 88 Her with M s = 5M , R s = 4R Q 
and T e = 12500K and with a disc defined by a = 2.5 
and (3 — 7.5. The disc has a free inner boundary and 
extends to r ~ 19. Plotted are the density profile E(r) of 
the unperturbed disc with the eigenfunctions E' and v' r 
for the m=l mode with P ~ 3.3 yr. The vertical scale is 
arbitrary. 

suiting disc profiles and calculated m=l eigen modes are 
plotted in Fig. 0] and respectively. The periods of the 
calculated prograde m= 1 modes are shown in Fig [S] For 
the disc model that touches the stellar surface ((3 = 0) 
we have calculated the m=l modes for the standard disc 
temperature Td = 2/3 T e and also for a cooler disc with 
Td = 6000 K. The disc model with free inner boundary 
{[3 = 7.5) exhibits prograde m=l modes with much lower 
oscillation periods and for significantly smaller values of 
the quadrupole factor &2 f 2 as compared to the discs with 
a rigid inner boundary. To compare these results with 
those obtained using the semi-analytic theory, we use the 
values of /C2/ 2 GM s /(c 2 R s ) above which prograde modes 
exist given in tabled To obtain values of k 2 f 2 we take 
GM S / (R s c 2 ) = 2 x 10 3 which is appropriate to the model 
with (a = 2.5, (3 = 7.5) and GM s /{R s c 2 ) = 2.8 x 10 3 for 
the model with (a = 2.5, ,9 = 0) and T d = 6000 K. These 
values correspond to the asymptotic limits for P — > 00 in 
Fig. H3 and are found to be 0.002 and 0.006, respectively. 
These results are consistent with the data shown in Fig. EI 
In particular they confirm the much larger critical value 
of fe/ 2 that applies in the case when a rigid wall inner 
boundary condition was applied to a distributed disc sur- 
face density profile. 
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Fig. 5. A model of 88 Her with a disc defined by a — 2.5 
and = with Td = 8333 K. The disc touches the stellar 
surface at r=l where it feels a rigid inner boundary (v' r = 
0), the free outer disc boundary is at r ~ 40. Plotted are 
the density profile E(r) of the unperturbed disc with the 
eigenfunctions £' and v' r for the m=l mode with P ~ 15 
yr. The vertical scale is arbitrary. 
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Fig. 6. The oscillation (or V/R) periods in years versus 
k.2p for the calculated (prograde) m=l modes for a model 
of 88 Her, as indicated in the figure. The results for the 
two disc models (for two different disc temperatures) with 
(3 — correspond to m = 1 modes with rigid wall inner 
boundary conditions. 

4. Conclusions 

We have investigated the effect of the rotational flattening 
of the rapidly spinning Be stars on low-frequency global 
m=l oscillations in their circumstcllar discs using both 
semi-analytic and numerical approaches. 

We considered both a hot, early type Be star sys- 
tem like 7 Cas and a cool late type system like 88 Her 
and found that if it is assumed there exists a (narrow) 
low surface density region or gap between the star and 
the disc, such that the disc inner boundary is free, low- 
frequency global m=l oscillations can be supported with 



periods in the observed range of the long-term V/R vari- 
ations of the Balmer emission lines for plausible values of 
the rotational factor ki] 2 (< 0.008). T herefore it is not 
necessary to introduce radiation forces (Okazaki 1997) to 
obtain prograde m=l modes that are contained in a oscil- 
lation cavity within a few stellar radii from a hot Be star 
surface. 

Our calculations suggest that the more compact cir- 
cumstcllar discs exhibit m=l modes with smaller oscilla- 
tion periods than the more extended discs for which the 
density maximum occurs slightly further out in the disc 
and for which the density decreases more slowly with in- 
creasing distance from the star. 

Of course in its present linear form our simple dynam- 
ical model cannot explain the full complexity of the ob- 
served variations in Be st ars, in particular non-l inear ef- 
fects need to be included ijFift & HarmaneJl2006l) espe - 
cially in enigmatic systems like 7 Cas ( Motch et al. 2005) , 
but it appears to lead to a natural explanation for the long 
term V/R variations in Be star systems, given the obser- 
vational fact that there are circumstcllar discs with slowly 
revolving m— 1 patterns. The model naturally acc ounts for 
the observed l|Berio et al.ll999UVakili et al.ll998h confine- 
ment and the prograde character of these patterns. How 
the circumstellar discs are formed and how the low fre- 
quency one armed density waves therein are excited still 
remains a puzzle. 
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